Dynamical surface structures in multi-particle-correlated surface growths 
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We investigate the scaling properties of the interface fluctuation width for the Q-mer and Q- 
particle-correlated deposition-evaporation models. These models are constrained with a global con- 
servation law that the particle number at each height is conserved modulo Q. In equilibrium, the 
stationary roughness is anomalous but universal with roughness exponent a — 1/3, while the early 
time evolution shows nonuniversal behavior with growth exponent j3 varying with models and Q. 
Nonequilibrium surfaces display diverse growing/stationary behavior. The Q-mer model shows a 
faceted structure, while the Q-particle-correlated model a macroscopically grooved structure. 
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I. INTRODUCTION 

Since the dynamical scaling theory for kinetic surface 
roughening |jj was suggested, numerous surface growth 
models have been studied M. The interface fluctuation 
width W developed during time t, starting from the flat 
surface, follows the scaling ansatz 



W ^ = La f^ 



(1) 



where L is linear size of the substrate. The scaling func- 
tion f(x) — > const, for x 3> 1 and f(x) ~ x 13 (j3 = a/zw) 
for x <C 1 pL^|- One-dimensional (ID) models with 
conventional local growth dynamics display the univer- 
sal roughness exponent a = 1/2, as exemplified in the 
Kadar-Parisi- Zhang (KPZ) universality class || and the 
Edwards- Wilkinson (EW) class {§. 

A ID surface can be mapped on the time trajectory of 
a particle in ID by identifying the height h(x) at each col- 
umn x with the particle position n(t) at time t — x. The 
step heights are bounded and uncorrelated beyond a fi- 
nite distance in the substrate direction. This implies that 
the particle performs a random walk without long-range 
diffusion and long-range temporal correlations, which is 
the characteristic of the normal random walks. There- 
fore the interface width scales with exponent a = 1/2, 
corresponding to the dispersion of normal random walks: 
An(t) -i 1 / 2 . 

Recently, Noh, Park, and den Nijs JB| introduced a 
dissociative dimer deposition-evaporation model with re- 
stricted solid-on-solid (RSOS) constraint. The surface 
grows and erodes by the deposition and evaporation of 
dimers aligned with the surface. These dimers dissoci- 
ate on the surface such that evaporating dimers do not 
necessarily consist of original partners. They found that 
the equilibrium surface is anomalously rough with a — 
0.29(4) and conjectured the exact value of a = 1/3 by 
exploring a close relationship between the dimer model 
and the even- visiting random walk model j^,^] . They also 
argued that the anomalous value of a = 1/3 is universal 
in generalized Q-mer type equilibrium models (Q > 2). 



The Q-mer aspect requires that the number of parti- 
cles at each surface height level must be conserved mod- 
ulo Q. The dissociative nature translates it into a non- 
local global constraint. So ID surfaces of Q-mer mod- 
els correspond to time trajectories of a random walker 
with a global constraint that the walker should visit each 
site (height level) Q multiple times before it terminates. 
The dispersion of Q-visiting random walkers (QVRW) 
is much more suppressed than the normal ones and be- 
haves asymptotically as An{t) ~ t 1 / 3 for Q > 2 §-§. 
This leads to the conjecture that the roughness exponent 
is a = 1/3 for general Q-mer type equilibrium models 
except for the monomer limit (a — 1/2 for Q = 1). 

However, the QVRW ensemble is not equivalent to 
the ensemble generated by the Q-mer dynamics. In 
fact, there exist infinitely many constants of motions 
in the Q-mer dynamics. For example, the Q-mer dy- 
namics conserves a local-type quantity. Upon depo- 
sition/evaporation of a Q-mer, the surface heights at 
Q consecutive columns change by one unit simultane- 
ously. It leads to conservation of k = 2n/Q compo- 
nent of the Fourier transformed surface height, defined 
as h(k) = Yl x =i ex P I~ ikx]h(x). Surface configurations 
with different values oih(k = 2n/Q) are dynamically dis- 
connected by the Q-mer dynamics. The phase space for 
the QVRW ensemble is divided into dynamically discon- 
nected sectors, labeled by h (k = 2n/Q). However, there 
is much more complex structure inside each sector. Each 
sector decomposes into infinitely many disconnected sub- 
sectors, which can be characterized by a nonlocal con- 
struct, i.e. so called irreducible string gj. Therefore, the 
Q-mer dynamics is strongly non-ergodic and the scaling 
property may be sector-dependent. 

Starting from a flat surface of size L as multiples of 
Q, one can explore only a subsector (null string) with 
h(k = 2tt/Q) = 0. Numerical results for dimer dynamics 
pj show strong finite size corrections and the extrapo- 
lation to infinite size limit overshoots significantly the 
expected value of a = 1/3 (5|, even though statistical er- 
rors are too large to conclude whether this overshooting 
is real. 
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This may be due to strong non-ergodic behavior of the 
Q-mer dynamics. In Ref. ||, Noh et al allow diffusion of a 
single particle on a terrace, which restores broken ergod- 
icity partially. Jumps across steps are forbidden. A cer- 
tain combination of deposition/evaporation and diffusion 
may be equivalent to a process of deposition/evaporation 
of Q particles not necessarily at Q consecutive columns 
but randomly placed on a terrace. It clearly does not 
conserve h(k = 2-k/Q) and opens some dynamic connec- 
tion links between sectors. Numerical results for Q = 2 
still show strong finite size corrections but the overshoot- 
ing is weakened and the estimated asymptotic value of 
a = 0.31(3) is a little bit closer to the conjectured value 
of 1/3 |. 

In this paper, we introduce a new model that fully re- 
stores the ergodicity. We allow deposition/evaporations 
of Q particles at randomly chosen Q columns with equal 
surface heights. These columns do not necessarily share 
the same terrace. These processes still conserve the 
global mod Q conservation at each height level and 
open all possible dynamic connections between sectors. 
We call this model as Q-particle-correlated deposition- 
evaporation model. The ensemble generated by this new 
type dynamics should be equivalent to the QVRW en- 
semble. As we shall see later, we find that the finite size 
corrections for Q = 2 are still considerable but the over- 
shooting is fairly reduced. Our estimate for a is now 
a ~ 0.33(3), which is fully consistent with the QVRW 
conjecture. 

We also investigated the scaling properties for Q-mer 
models and Q-particle correlated models for Q = 3 and 
4. Interestingly, the asymptotic scaling seems to set in 
much earlier than Q = 2 models and the estimated val- 
ues of a are very close to 1/3 in both models. These 
results may suggest that the phase space structures for 
Q = 2 and Q > 3 bear an intrinsic difference like in the 
monolayer version of the Q-mer models [|| . 

Time-dependent behavior in early time regime is much 
more complex and sector-dependent. Various interest- 
ing dynamic properties have been reported by Grynberg 
in the body-centered solid-on-solid (BCSOS) version of 
Q-mer models [ ff0f . In this paper, we focus only on a 
subsector (null string) containing a flat surface. Our nu- 
merical results show that the scaling exponent j3 in Eq. (jl]) 
varies slightly but definitely with Q in both Q-mer and Q- 
particle correlated models. Our estimated values of (3 for 
the Q-mer models are consistent with those for the BC- 
SOS version, but different from those for the Q-partAcile 
correlated models. It may reflect non-ergodic nature of 
the Q-mer dynamics. In contrast to the universal scaling 
property in the steady state (universal a), the dynamic 
properties characterized by (3 vary with models and Q. 

Excursion off equilibrium immediately changes the sur- 
face structure drastically ||. Without balance between 
deposition and evaporation, the surface grows or erodes 
indefinitely with time. In the Q-mev models, these 



nonequilibrium surfaces turn out to be always faceted 
due to spontaneous formation of pinning valleys or pin- 
ning hill tops ||. In finite systems, the surface displays 
a repeated jerky motion in time, i.e., alternating series 
of locking in a facet shape and unlocking/evolving into 
another facet shape. Time scale of this jerky motion is 
dominated by duration of a facet shape, which diverges 
exponentially with system size. The global facet struc- 
ture in the steady state assures that the interface width 
W scales as system size L (a = 1). Along with the ex- 
ponential time scale, one can predict that W ~ logt for 
early time regime. These results are confirmed numeri- 
cally. 

In case of the Q-particle-correlated models, nonequi- 
librium surfaces are not sharply faceted but display a 
global grooved structure with considerable local fluctua- 
tions. This surface is similar to those found in the con- 
served RSOS (CRSOS) model Q. The roughness expo- 
nent is found to be a = 1.00(5) for all Q > 2. Early time 
growth behavior is governed by a power law with nonuni- 
versal exponent /3 varying with Q and the magnitude of 
deposition/evaporation bias. 

Outline of this paper is as follows. In Sec. II, we in- 
troduce the Q-particle-correlated models as well as the 
Q-mer models. In Sec. Ill and IV, we present numeri- 
cal results for the scaling properties of the equilibrium 
and nonequilibrium surfaces, respectively. Similarities 
and differences between two models are discussed. We 
conclude with a brief summary and discussion in Sec. V. 



II. MODELS 

Consider the ID surface configurations described by 
integer height variables h(x) = 0, ±1, ±2, ... with x — 
1,...,L. They are subject to the RSOS constraint, 
h(x) — h(x + 1) = 0, ±1, and periodic boundary con- 
ditions, h(x + L) = h(x). 

The growth rule for the Q-particle-correlated 
deposition-evaporation model is as follows. First, select 
Q columns {x±, x%, xq} randomly. These columns 
do not have to be adjacent to each other. If the 
heights {h(xj)} of the selected columns are all equal, 
i.e., h(xi) = ... = h(xq), then simultaneous deposition 
of Q particles, h{xi) — > h{xi) + 1 for i = 1,..,Q, is 
attempted with probability p, or simultaneous evapora- 
tion, h(xi) — ► h{xi) — 1 for i = 1, with probability 
q = 1 — p. Any attempt is rejected if it would result in 
violating the RSOS constraint. If any selected column 
height differs from other selected one, no dynamics oc- 
curs and we select a new set of Q columns randomly. 
After L such selection processes, Monte Carlo time t is 
incremented by one unit. 

The equal height condition, h(x\) = ... = h(xQ), con- 
strains a set of selected columns and generates correlated 
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growth. Without this condition, any set of randomly se- 
lected columns can evolve and the ordinary monomer- 
type RSOS model is recovered. One can consider a more 
restrictive constraint on selected columns. An interest- 
ing and limiting case is to require the immediate adja- 
cency of selected columns. This is the dissociative Q-mcr 
deposition-evaporation model. In this Q-mer model, we 
attempt to select Q consecutive columns only. The other 
evolution rules are equivalent to the Q-particle-correlatcd 
model. The generalization of both models to higher di- 
mensions should be straightforward and will be discussed 
elsewhere Q. In this paper, we focus on ID cases only. 



III. EQUILIBRIUM SURFACES 

With balance between deposition and evaporation pro- 
cesses (p — q), the detailed balance condition is satisfied. 
In the stationary state, the Gibbs-type equilibrium sur- 
face is obtained. In this section, we investigate the scal- 
ing property of the equilibrium surfaces for the Q-mer 
and Q-particle-correlated models and also their dynamic 
property approaching to equilibrium. We measure the 
interface fluctuation width W as 



W\L,t) = W' 



h(x, i) 



1 L 



(2) 



where L is the substrate size. We start with a flat sur- 
face of size L as a multiple of Q, otherwise the surface is 
forcefully pinned to the initial position. 



A. stationary state roughness 

For finite size scaling analysis, we run Monte Carlo 
simulations for system size L = 2™ (n = 4, 10) in case 
of Q = 2 and 4, and L = 2 x 3™ (n = 2, ...,6) in case 
of Q = 3. The interface width is measured and averaged 
over at least 100 independent Monte Carlo runs for each 
system size. It grows in early time regime and saturates 
to a finite value, which depends on size L. We again 
average over data in the saturated regime (t 3> L zw ) to 
estimate the stationary value of the width, W(L, t = oo). 

As explained in Eq. (|l|), W(L,t = oo) scales as L a in 
the large L limit. For efficient data analysis, we introduce 
effective exponents a e f / as function of L 



aeff(L) 



In W(2L, t = oo) - In W(L, t = oo) 
ln(2L) -InL ' 



(3) 



The value of a can be retrieved by taking the L — > oo 
limit for a e ff(L). 

Effective exponents for the Q-mer model are plotted 
against 1/L in Fig. 1. For the dimer case, one can see 
usual 1/L type finite size corrections for small L, but 



strong crossover corrections near L = 2 8 . The extrapola- 
tion of a e f f to the infinite size limit seems to overshoot 
significantly the expected value of a = 1/3. However, 
the statistical errors are too large to conclude that our 
numerical results exclude the possibility of a — 1/3. In 
fact, it is quite difficult to reduce statistical errors con- 
siderably with moderate computing usage, because W is 
very small even for large L (less than 4 for L = 2 10 ). 
We estimate a = 0.29(4) for the dimer model. For Q = 3 
and 4, we do not find strong crossover corrections and our 
numerical data converge to a = 1/3 very nicely. We es- 
timate a = 0.32(4) for the trimer model and a = 0.33(3) 
for the 4-mer model. 

Data for the Q-particle-correlated models are shown 
in Fig. 2. For Q = 2, we find again strong crossover 
corrections, but the overshooting is fairly reduced. We 
estimate a = 0.32(3), 0.33(2) and 0.32(4) for Q = 2,3, 
and 4, respectively. These results are fully consistent 
with the QVRW conjecture. 



B. Time dependent behavior 

In this subsection, we investigate the early time be- 
havior of the interface width. We measure the interface 
width for system size L = 10 4 for Q = 2 and 4, and 
L = 2 x 3 6 for Q = 3 up to t = 10 5 - 5 x 10 6 and aver- 
age over at least 100 independent samples. In early time 
regime (t <C L zw ), it grows algebraically as W ~ t 13 , see 
Eq. ([!]). To estimate the value of (3, we introduce effective 
exponents (3 e f t as a function of t 



lnW(t) -1nW(t/10) 
lni - ln(t/10) ' 



(4) 



The value of (3 can be retrieved from the long time limit 
for (3 e ff within t <C L Zw , where finite size effects are 
negligible. 

Effective exponents for the Q-mer model and the Q- 
particle-correlated model are plotted against t in Fig. 3 
and 4, respectively. After initial transients, f3 e ff become 
stabilized around an asymptotic value. We estimate (3 
by averaging f3 e ff in the stabilized regime. As seen in 
Fig. 3, the estimated values for (3 are 0.108(1) for the 
dimer, 0.100(1) for the trimer, and 0.098(3) for the 4- 
mer model. Surprisingly, the value of (3 decreases slightly 
but definitely with Q. Our finding is consistent with the 
result for the BCSOS type Q-mer model [@. 

The Q-particle correlated models show even strong 
variance of (3 on Q. ^From Fig. 4, we estimate (3 = 
0.194(1), 0.140(2), 0.097(2) for Q = 2, 3, and 4. Our re- 
sults suggest that (3 continuously varies and approaches 
zero in the large Q limit. 

The steady state scaling is quite universal and does 
not depend on Q{> 2). This universality has been well 
established through a mapping to the Q-visiting random 
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walks Q. In contrast, our numerical results strongly sup- 
port nonuniversal dynamic scaling depending on Q. As 
yet, we do not have any reasonable explanation for this 
type of nonstandard scaling behavior. 

It is not quite surprising that the Q-mer and Q- 
particle-correlated models display different dynamic scal- 
ing behaviors. The Q-mcr dynamics are not ergodic, in 
contrast to the Q-particle-correlated dynamics. The non- 
ergodicity slows down the dynamics and the growth ex- 
ponent should be smaller, which is consistent with our 
numerical finding. 

It is interesting to see that initial transients for the 
Q-mer and the Q-particle-correlated models are com- 
pletely upside down. We can give a reasonable explana- 
tion for this. In case of the Q-mer equilibrium model, it is 
clear that the early time behavior resembles the Edward- 
Wilkinson type (monomer model), because the system 
does not feel the Q-mer aspect initially before pinning 
valleys/hilltops start to form. So, /3 e // should decrease 
from Pew = 0.25 to converge to the asymptotic value. 
The Q-mer aspect should appear earlier for larger Q with 
the same substrate size. We expect that the initial tran- 
sients survive for shorter period for larger Q, which can 
be seen in Fig. 3. 

In case of the Q-particle-correlated model, it is very 
difficult to grow the surface at the beginning. As Q par- 
ticles deposit on randomly placed Q sites, it is difficult to 
form aggregated islands on which one can deposit more 
particles to grow the surface. So in the very early time 
regime, the surface would grow in a layer-by-layer fash- 
ion. So P e ff will start from zero and increase to the 
asymptotic value. Here, the initial transients last longer 
for larger Q. The rejection rate for the Q particle deposi- 
tion/evaporation on the randomly selected sites increases 
exponentially with Q. In fact, the initial transient be- 
havior is so huge for Q = 4 that we use a rather small 
substrate size L = 10 3 in Fig. 4. 

IV. NONEQUILIBRIUM SURFACES 

Without balance between deposition and evaporation 
(p 7^ <?)j the surface may grow (p > 1/2) or erode 
{p < 1/2) indefinitely with time, if one starts with a flat 
surface of size I as a multiple of Q. Due to the time re- 
versal symmetry, the growing surface at deposition prob- 
ability p should be identical to the inverted eroding sur- 
face at 1 — p. In the monomer version (Q = 1), it is 
well known that the unbalance (p ^ q) becomes a rele- 
vant perturbation to the EW fixed point and drives the 
system into the KPZ universality class. In this section, 
we investigate the scaling property of the noncquilibrium 
growing/eroding surfaces for Q > 2. 

The nonequilibrium surfaces for the Q-mer model and 
the Q-particle-correlated model show completely distinct 
characteristics in both dynamic evolution and stationary 



morphology. The former ones always facet, while the lat- 
ter ones form global grooved structures. For clarity, we 
discuss these two models separately in this section. 



A. Q-mer model 

The noncquilibrium surfaces of the Q-mer model dis- 
play sharply faceted structures. Fig. 5 shows the time 
evolution of surfaces for the dimer model at p — 0.6 and 
0.1 in a typical simulation sample. In the growing sur- 
faces, the faceting is caused by the spontaneous forma- 
tion of pinning valleys. For the dimer model, any flat 
segment with odd size acts as the nucleus of a pinning 
valley. Such valleys can not be filled up by deposition 
of dimers and one empty blocked site always remain at 
the valley bottom. Pinning valleys can annihilate only 
in pairs, and only if their valley bottoms are at the same 
height. Therefore, the only way to remove pinning val- 
leys and grow the surface is the evaporation of the entire 
hill between two pinning valleys. 

In finite systems, the surface moves like shock waves. 
The initially flat surface evolves into a faceted shape 
with only two remaining pinning valleys. The valleys 
are sharp and the hilltops are rounded with a typical size 
£o- Approaching the equilibrium point, £o diverges and 
the facets disappear. The annihilation time of this last 
pair scales exponentially with system size. After its an- 
nihilation, the surface grows into another faceted shape 
very fast and the whole evolution repeats itself. Details 
of this motion have been reported in Ref . || . 

For general Q-mer models, there exist Q — 1 different 
kinds of pinning valleys with n empty blocked sites at the 
valley bottom (n — 1,2, ...,Q — 1). Two pinning valleys 
with n and m empty blocked sites merge into a pinning 
valley with n + m empty blocked sites. When n + m is a 
multiple of Q, the valley can be filled up and the surface 
can grow. Even though the surface evolution details are 
different for each Q, the basic evolution characteristics 
are equivalent, i.e., a repeated jerky motion of locking 
into a facet shape and sudden evolving into another facet 
shape. Time scale of this jerky motion is dominated by 
exponentially long duration of a facet shape. If one starts 
with a flat surface of size L that is not a multiple of Q, at 
least one pinning valley always survives and the surface 
is pinned to the initial position. 

The interface width W of the faceted surface in Fig. 5 
should be proportional to the substrate size L. There- 
fore, the stationary roughness exponent is trivially a = 1. 
Effective exponents ct e ff for the Q-mer model are plot- 
ted against 1/L in Fig. 6. As expected, <x e ff converge to 
1 for all Q. For small sizes, finite size effects are present 
due to finite duration time of the facet shape. 

The early time growing behavior also can be predicted 
from annealing dynamics of pinning valleys. Starting 
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from a flat surface, the surface grows fast at the begin- 
ning and many pinning valleys appear. The presence of 
pinning valleys slows down the growth and the surface 
may grow only when two pinning valleys merge together 
by annealing out the hill between them. As seen in Fig. 5, 
only a small region on the hilltops can be active by depo- 
sition/evaporation processes. Near the valley, the ramps 
are inactive due to the RSOS constraint. So, the anneal- 
ing time of the hill should scale exponentially with lateral 
size of the hill. After this time scale, the surface advances 
fast locally up to the order of the hill height. In the 
faceted structure, the hill height should be proportional 
to the hill size. In this point of view, one can conclude 
that the local surface fluctuation should grow logarith- 
mically with time. As any other time scale is present, we 
expect that W ~ logi in the early time regime |j 10 1. 

In Fig. 7, we plot W against logt. For all Q, it 
shows a nice linear behavior. For clarity, we assume that 
W ~ [logi] x and plot effective exponents Xefft defined 
similarly in Eq. (|j), versus t in the inset. They con- 
verge rather slowly but nicely to 1. We estimate that 
X = 1.1(1). 



B. Q-particle-correlated model 

Surface morphology for the Q-particle-correlatcd 
model is very surprising. Fig. 8 shows the time evolu- 
tion of surfaces for Q = 2 at p = 0.6 and 0.1 in a typical 
simulation sample. In contrast to the Q-mer model, the 
surfaces are not faceted but show a global grooved struc- 
ture with considerable local fluctuations. 

Starting from a flat surface, small scale hills and valleys 
start to form and merge into a few macroscopic hills and 
valleys. These macroscopic hills/valleys contain consid- 
erable roughness in microscopic scale. Especially the hill- 
tops of the growing surface are broad and rough with high 
activity of deposition/evaporation processes. In contrast, 
the valley bottoms are rather localized and their ramps 
are quite steep without much activity. However, they 
are not extremely localized or sharp as in the faceted 
structure of the Q-mer model. The valley bottoms do 
not last long exponentially with system size, but only al- 
gebraically. Therefore, there is no well-defined repeated 
motion of the growing surface. A few macroscopic valleys 
merge and split stochastically in the steady state regime. 

In Fig. 9, we plot a e ff against 1/L for Q — 2,3, and 
4. We estimate that a = 1.00(5) for all Q. In contrast 
to the Q-mer model, the early time growth behavior for 
the Q-particle-correlated model is governed by a power 
law with exponent /3 varying with Q and p. We plot (3 e f f 
against t at p = 0.6 in Fig. 10. Our estimates at p = 0.6 
are (5 = 0.460(1), 0.325(1), and 0.214(2) for Q = 2, 3 and 
4. The value of (3 increases with increasing p (or depo- 
sition/evaporation bias). The power-law type growth is 
consistent with our groove formation picture. The ramps 



near the valley are still active, so the annealing time of 
the hill between two valleys does not scale exponentially 
but algebraically. 

A similar grooved structure could be found in the so 
called conserved RSOS (CRSOS) model O and other 
models with roughness exponent a > 1 |2]7 It may be 
quite interesting to investigate similarities and differences 
among these models with grooved structures. For exam- 
ple, the grooved structure in the CRSOS model displays a 
sharp peak and a rounded valley, which is just the reverse 
of that in our model. Moreover, the physical origin for 
occurrence of macroscopic grooves is not well explored. 
More detailed investigation is left for the future study. 



V. SUMMARY AND DISCUSSION 

We investigated the scaling properties of the inter- 
face fluctuation width for the Q-mer and Q-particle- 
corrclatcd deposition-evaporation models. Both models 
are constrained with a global conservation law that the 
particle number is conserved modulo Q at each height 
level. A ID surface of these models can be mapped on a 
random walk trajectory with Q-visiting global constraint. 

When deposition and evaporation are balanced, one 
can obtain Gibbs-type equilibrium surfaces in the sta- 
tionary state. Then, the dispersion of the Q-visiting 
random walks (QVRW) can translate into the station- 
ary roughness of equilibrium surfaces. However, the Q- 
mer dynamics is not ergodic and can not generate all 
possible QVRW configurations. This non-ergodicity may 
generate huge corrections to scaling and hinder the nu- 
merical analysis for the stationary roughness. Especially 
for the dimer model, we found a significant overshooting 
of roughness exponent a in finite-size-scaling analysis. In 
this paper, we introduce the Q-particle-correlated model 
that fully restores the ergodicity. Our numerical study 
showed that the overshooting is fairly reduced and our 
estimate for a is fully consistent with the QVRW result. 

We also studied the time-dependent behavior in the 
early time regime. Interestingly, the growth exponent (3 
varies with models and Q (see Table 1), in contrast to 
the robustness of a. Therefore, these models correspond 
to a series of new universality classes. 

When the evolution dynamics is biased, the surface can 
grow or erode indefinitely with time. As the detailed bal- 
ance is violated, the steady-state surface configurations 
do not form a Gibbs-type ensemble. We numerically in- 
vestigated the scaling properties of these nonequilibrium 
surfaces. The Q-mer models display shock- wavelike evo- 
lutions and finally form a macroscopic facet, while the Q- 
particle-correlated models follow an ordinary power-law 
type growth and form a macroscopic grooved structure. 
Like in the equilibrium cases, the steady state rough- 
ness do not depend on Q, but the growth exponent for 
the Q-particle-corrclatcd models varies with Q and the 
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deposition/evaporation bias. These values are listed in 
Table 1. 

In summary, we found a series of new universality 
classes in deposition-evaporation models with a global 
conservation law. Even though the stationary roughness 
in equilibrium can be explained in terms of QVRW mod- 
els, the growth diversity is not fully understood. The 
scaling properties for growing/eroding surfaces for the 
Q-particle-correlated models are surprising and wait for 
a reasonable intuitive explanation. Generalization to 
higher dimensions are under current investigation. 
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TABLE I. Roughness exponent a and growth exponent /3 for the Q-mer and Q-particle-correlated (Q-PC) models. For 
nonequilibrium Q-PC models, /3 also varies with p. 
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Q = 2 


= 3 


Q = 4 


Q = 2 


Q = 3 


= 4 
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0.29(4) 


0.32(4) 


0.33(3) 


0.108(1) 


0.100(1) 


0.098(3) 


(P = <?) 


Q-PC 


0.32(3) 


0.33(2) 


0.32(4) 
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FIG. 1. Effective stationary roughness exponents a e ff versus 1/L for the Q-mer deposition-evaporation models at 
p = l/2(equilibrium surfaces). The horizontal broken lines represent a e // = 1/3. 



7 



0.5 



0.4 



0.3 




0.00 



0.02 



0.04 



0.06 



1 /3 



1 I L 



FIG. 2. Effective stationary roughness exponents a e ff versus 1/L for the Q-particle-correlated deposition-evaporation 
models at p = 1/2. The horizontal broken lines represent a e ff = 1/3. 
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FIG. 3. Effective growing exponents f3 e ff versus time t for the Q-mer models at p — 1/2. The horizontal lines are 
(3 eff = 0.108, 0.100, and 0.098 for Q = 2,3, and 4, respectively. 
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FIG. 4. Effective growing exponents /3 e // versus time t for the Q-particle-correlated models at p — 1/2. The horizontal lines 
are /3 e ff = 0.194, 0.140, and 0.097 for Q — 2, 3, and 4, respectively. We take a rather small substrate size L — 10 3 for Q = 4, 
due to huge initial transient behavior. 
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FIG. 5. Time evolutions of surfaces for dimer growth model at p — 0.6 and p = 0.1 in a typical simulation sample. 
The numbers in the bottom of figures denote Monte Carlo times when each surface configuration is taken. Typical surface 
configuration with maximum interface width are drawn in the right for Q = 2, 3, and 4. 
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FIG. 6. Effective stationary roughness exponents ot e ff versus 1/L for the Q-mer models at p = 0.6. The horizontal broken 
lines represent a e // = 1. 
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M 

FIG. 7. Interface width W against \ogt for the Q-mer models at p = 0.6 in t <S L Zw . Assuming that W ~ [logt] x , we plot 
effective exponents Xeff versus t in the inset. Xeff converge rather slowly but nicely to 1. 
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FIG. 9. Effective stationary roughness exponents a e // versus 1/L for the Q-particle-correlated models at p — 0.6. The 
horizontal broken lines represent a e // = 1- 
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